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ABSTRACT 

Apparent redshift dependence of the jet opening angles (#j) of gamma-ray 
bursts (GRBs) is observed from current GRB sample. We investigate whether 
this dependence can be explained with instrumental selection effects and obser- 
vational biases by a bootstrapping method. Assuming that (1) the GRB rate 
follows the star formation history and the cosmic metallicity history and (2) the 
intrinsic distributions of the jet-corrected luminosity (L 7 ) and 9- } are a Gaussian 
or a power-law function, we generate a mock Swift/BAT sample by considering 
various instrumental selection effects, including the flux threshold and the trigger 
probability of BAT, the probabilities of a GRB jet pointing to the instrument 
solid angle and the probability of redshift measurement. Our results well repro- 
duce the observed 6- } — z dependence. We find that in case of L 7 oc # 2 good 
consistency between the mock and observed samples can be obtained, indicating 
that both L 7 and 9- } are degenerate for a flux-limited sample. The parameter set 
(L 7 , 9j) = (4.9 x 10 49 erg s _1 , 0.054rad) gives the best consistency for the current 
Swift GRB sample. Considering the beaming effect, the derived intrinsic local 
GRB rate accordingly is 2.85 x 10 2 Gpc -3 yr -1 , inferring that ~ 0.59% of Type 
Ib/c SNe may be accompanied by a GRB. 

Subject headings: Gamma-ray burst: general-Methods: statistical-Stars: forma- 
tion 



1. Introduction 

Gamma-ray bursts (GRBs) are the most luminous explosions in the Universe. They are 
detected from deep to our local Universe (Zhang & Meszaros 2004). The measured redshifts 
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of current GRB sample range from 0.0085 (GRB 980425; Tinney et al. 1998) to 8.2 (GRB 
090423; Tanvir et al. 2009) 1 . It is generally believed that the births of long GRBs follow 
the star formation history of the universe (e.g., Totani 1997; Paczynski 1998; Bromm & 
Loeb 2002; Lin et al. 2004), and they may be promising probes for cosmology and galaxy 
evolution (e.g., Dai et al. 2004; Ghirlanda et al. 2004; Liang et al. 2005; Yiiksel & Kistler 
2007; Kistler et al. 2008, 2009; Salvaterra et al. 2008; Wang & Dai 2009). However, it is 
still unclear if GRBs experience any sort of cosmic evolution. 

With high sensitivity of the burst alert telescope (BAT) and the promptly slewing 
capability of the narrow-field X-ray telescope (XRT) onboard the Swift satellite (Gehrels et 
al. 2004), the number of measured high-z GRBs has increased rapidly over the last 6.5 
years. It is found that long GRBs do not follow star formation rate (SFR) unbiasedly (e.g., 
Kistler et al. 2008). The observed GRB rate is higher than the SFR at z ~ 4, which may 
be due to that the GRB rate is also related to the comic metallicity history (Kistler et al. 
2008; Li 2008; Qin et al. 2010; Virgili et al. 2011) or due to a higher efficiency of long 
GRB production rate by massive stars at high- z (Daigne et al. 2006). It was also suggested 
that the GRB luminosity may evolve with redshift, being more luminous in the past, which 
is simplified as L oc L (l + z) 6 , where L is the luminosity of local GRBs (Lloyd- Ronning, 
Fryer & Ramirez-Ruiz 2002; Firmani et al. 2004; Kocevski & Liang 2006; Salvaterra et al. 
2008, 2009). 

It is believed that GRBs should be collimated into narrow jets (Rhoads 1997). The 
collimation of GRB jet is critical to understand the physics of GRBs, such as their central 
engine, intrinsic GRB rate, and true gamma-ray energy explosion (e.g., Zhang et al. 2004). 
The observed steepening with a decay slope of ~ 2 in the late afterglow lightcurves is possibly 
due to the jet effect, in which a break would happen when the bulk Lorentz factor of the 
outflow (r) slows down to T < Qj 1 , where 9 } is the half opening angle of the jet (Rhoads 
1997; Sari et al. 1999). With the observed jet break time tj, the 9- } can be calculated (Rhoads 
1997; Chevalier & Li 2000). Frail et al. (2001) derived the geometrically-corrected gamma- 
ray energy (i? 7 ) of a small GRB sample and proposed an idea of a standard energy reservoir 
in GRB jets based on the narrow distribution of E 1 (see also Berger et al. 2003 for the 
same conclusion for the jet-corrected luminosity (L 7 ). However, this result seems to be a 
sample selection effect. The E 1 distribution is significantly broadened with the current GRB 
sample, although it is still log-normal (Liang et al. 2008; Racusin et al. 2010). It was also 
proposed that GRB jet opening angles may experience similar kind of evolution as GRB 
luminosity (Wei et al. 2003). Estimated the 9 } of the BATSE GRB sample with an empirical 



1 A photometric redshift of 9.4 for GRB 090429B was reported by Cucchiara et al. (2011) 
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relation based on the Yonetoku and Ghirlanda relations (Ghirlanda et al. 2004; Yonetuku 
et al. 2004), Yonetoku et al. (2005) argued that 0j evolves as 9j oc (1 + z)~ 0A5 . 

The evolution features of GRB luminosity and jet opening angle as well as the GRB 
jet energy distribution are based on the statistics for the current observational samples. 
However, as well known, observational samples highly suffers instrumental selection effects 
and observations biases (Bloom 2003). It is difficult to reveal the intrinsic distributions and 
evolution features of GRBs from these samples. For example, a GRB with a peak photon 
flux lower than an instrumental threshold cannot trigger the instrument. Even the flux is 
over the threshold, the trigger probability of a GRB strongly depends on its photon flux (Qin 
et al. 2010). The probability of a GRB jet pointing to the solid angle of an instrument is 
different for GRBs with different 9 y In addition, E 7 (or L 7 ) and 9j are not direct observable. 
The derivations of these parameters are mainly based on the observed flux, redshift, and t y 
However, both E y (or L 7 ) and 9 S are degenerate for a flux-limited sample. This means that 
an intrinsically dimmer GRB (with lower E 1 or L 7 ) may be detected by a given instrument 
in case of a smaller 9 y 

In this paper we investigate whether the observed redshift dependences of jet opening 
angles and isotropic gamma-ray luminosity of long GRBs can be explained with instrumental 
selection effects and observational biases. Based on our results we also estimate the local 
long GRB rate by considering the beaming effect. We search for GRBs whose jet break 
times in afterglow lightcurves are reported in literature. The data are reported in Section2. 
Our models and results for Swift/BAT sample are reported in section 3. Conclusions and 
Discussion are presented in Section 4. Throughout we use cosmological parameters Hq = 
71km s -1 Mpc"\ VL M = 0.3, fi A = 0.7. 



2. Data 

Our sample includes all GRBs whose jet break times (tj) are measured in the radio, 
optical, and X-ray afterglow lightcurves, regardless the t- s are achromatic, or just detected 
only in one band. Note that a tentative achromatic jet breaks in the radio, optical and X-ray 
afterglow lightcurves are available for three GRBs. We take the t- s of these GRBs measured 
in the optical afterglow lightcurves. Seventy-seven GRBs are included in our sample. They 
are summarized in Table 1. Since different cosmological parameters and medium density 
surrounding the bursts are used in calculations of 9- } in literature, we re-derive £i so , 9 y and 
E 7 in a consistent method with the same cosmological parameters and the same jet model 
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for a constant medium density (Rhoads 1997, Sari et al. 1999, Frail et al. 2001) 

S = .057 r ad (r i-P(i±f )- 3 / 8 ( I ^)-V 8( i )1 /s ( _^_ )1/8 , (1) 

where n is the ambient medium density in unit of cm -3 and r/ 7 is the efficiency of the fireball 
in converting the energy in the ejecta into gamma-rays. We take r? 7 = 0.2 and n ~ 0.1 cm -3 
in our calculations. The E iso is calculated with 



E^=^K (2) 

where S 1 is the observed gamma-ray fluence, D L is the luminosity distance at redshift z, 
and k is a factor to correct the observed gamma-ray energy in a given band pass to a broad 
band (i.e., 1 — 10 4 keV in this analysis) with the observed GRB spectra (Bloom et al. 2001). 
We collect the spectral parameters for the burst in our sample from literature. It is well 
known that the GRB spectrum is well fit with the Band function (Band et al. 1993). Since 
the narrowness of the Swift/BAT band, the spectra of most Swift GRBs in our sample are 
adequately fit with a single power-law, N oc E~ r (Zhang et al. 2007; Sakamoto et al. 
2009). If the parameters of the Band function are unavailable for the Swift GRBs, we use 
an empirical relation between the peak energy of the vf v spectrum and the photon index T, 
i.e. log E p = (2.76 ± 0.07) - (3.61 ± 0.26) logT (Zhang et al. 2007), to estimate the E p and 
take the low and high photon indices as a = — 1.1 and /3 = —2.2 (Preece et al. 2000; Kaneko 
et al. 2006). The geometrically-corrected gamma-ray energy of a GRB jet can be obtained 
with 

£ 7 = (l-cos0j)£ iso . (3) 

The derived Q- } and E 1 are reported in Table 1. It is found that their distributions are 
log-normal, as shown in Figure 1. The fits with a Gaussian function yield log(#j/rad) = 
(-1.31 ±0.24) (la) and log(£ 7 /erg) = (50.07±0.91) (la), respectively. The 9 i as a function 
of z is shown in Figure 2. A tentative anti-correlation is observed. The best fit gives 

log#j = (-0.90 ±0.09) - (0.94 ± 0.19) log(l + z), (4) 

with a Spearman correlation coefficient of 0.55, a standard deviation of 0.30, and a chance 
probability p < 10~ 4 for N = 77. 



3. Analysis of the Instrumental Selection Effects 

It is unclear whether the observed 9j — z dependence shown above is due to observational 
selection effects or due to an intrinsic cosmological evolution feature of GRBs. In this section, 
we make bootstrap analysis of the instrumental selection effects on this dependence. 
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3.1. Models 

The number density of GRBs at redshift z ~ z + dz is given by 

t t \ dN R GRB{z)dV{z) 
U ^ ^ = dzlL^T = 1 + z dz ^ (5) 

where Rgrb(z) is the co-moving GRB rate as a function of z, the factor (1 + z)~ l accounts 
for the cosmological time dilation of the observed rate, dT is the time interval, $(£ 7 ) is the 
luminosity function of GRB jets and dV(z)/dz is the co-moving volume element. We assume 
-Rgrb traces the star formation rate and metallicity history, reading as (Kistler et al. 2008; 
Li 2008; Qin et al. 2010; Virgili et al. 2011) 

Rgrb(z) = p o ^SFR(^)0(e, z), (6) 

where po is the local GRB rate, 0(e, z) is the fractional mass density belonging to metallicity 
below eZ & at a given z (Z Q is the solar metal abundance) and e is determined by the 
metallicity threshold for the production of GRBs. In our analysis, the Rsfr(z) is taken as 

Rsfr(z) oc <^ 44 , (7) 

^ ( 1 "i" ^peak ) i Z > 2p Ca k 

where z pea k = 1 (Rowan-Robinson 1999). The 9(e, z) is parameterized as (Hopkins & Beacom 
2006; Langer & Norman 2006) 

f(a + 2,e*10°- 15 *) 

0(e ' 2) = fvT2) ' (8) 

where a and (5 are the low and high spectral indices of the prompt gamma-rays, T(a,x) and 
T(x) are the incomplete and complete gamma function (Langer & Norman 2006; Kistler et 
al. 2008). We adopt the jet-corrected luminosity function $(£ 7 ) as a log-normal distribution 
function (Matsubayashi et al. 2005), 



*(L T ) = ^exp[- " ^-' 0gLc)2 l. (9) 



2X01, 1 2<rJ 



7 



where is a normalized constant, L c is the center value and is the width of the 
distribution. The intrinsic distribution of Q- } is unknown. We also take it as a log-normal 
distribution 2 , 

nOj) = tk= ex p[ -^—2 ]> ( 10 ) 



2 The case of a broken power-law and a power-law distributions for L 7 and 0j, respectively, is discussed 
in section 4. 
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where Aq is a normalized constant, C is the center value and oq is the width of the distribu- 
tion. The observed peak flux in a given instrumental band for a GRB that has a half-opening 
angle Oj and jet-corrected luminosity P 7 at redshift z is derived from 

P(z, L 7 , Oj) = 4nD 2 L ^ k = 4^1(^(1-008^)^' (U) 

where L = L 7 / (1 — cos 

The algorithm of the Swift/BAT trigger has two modes, the count rate trigger and image 
trigger (Fenimore et al. 2003; Sakamoto et al. 2008, 2011). Rate triggers are measured on 
different time scales (4 msec - 64 sec), with a single or several backgrounds. Image triggers are 
found by summing images over various timescales and searching for un-catalogued sources. 
Band (2006) compared the threshold of Swift/BAT in count rate mode with other GRB 
missions. The image trigger depends on the duration of a burst. Sakamoto et al. (2008) 
parameterized the threshold of BAT as a function of burst duration. With the trigged and off- 
line scanned samples of CGRO/BATSE, Qin et al. (2010) simplified the trigger probability 
as a function of the observed peak flux for BATSE. They took the same trigger probability 
for BAT to make bootstrap analysis for the BAT sample. In this analysis, we take a modified 
form of the trigger probability as proposed by Qin et al. (2010). As shown by Sakamoto et 
al. (2008), the trigger probability of the image trigger mode for those long GRBs with a low 
peak photon flux would be increased. In order to increase the trigger probability of long and 
weak GRBs that have a peak flux near the BAT threshold with the image trigger mode, we 
take the BAT trigger probability as a function of P as following 3 , 

, N [P 2 P < 0.45 ph/cm 2 /s , 

[0.67(1.0 - 0.40/P) - 52 , P>0.45 ph/cm 2 /s 

The probability of redshift measurement for a GRB is complicated and depends on many 
parameters (Fynbo et al. 2009). It is difficult to parameterize the redshift measurement. 
Qin et al. (2010) got a weak dependence of the probability on the observed peak flux , which 
is quoted as following, 

7] Z (P) = 0.26 + 0.032e L611ogP . (13) 

The probability of the alignment for a GRB with jet opening angle 0j to an instrument with 
a solid angle Q for the aperture flux is estimated with 

Va(0 j ) = ^(l-cos0 J ). (14) 



3 Notc that we focus on the redshift-known sample in this analysis. These GRBs are usually bright. 
Therefore, some uncertainties of the trigger probability of weak GRBs near the instrumental threshold 
would be not significantly affected our analysis results. 
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Therefore, the number of GRBs with redshift measurement that trigger an instrument 
in an operation period T can be calculated with 

max 

n = t ^{e j )de j / ^(p)^ 2 (p)^(^)$(L 7 )rfL 7 

J ^^7, min / 1 K\ 

x [ z — Rgrb{z) dV(z) dz 
Jo 1 + z dz 

where z max is determined by the instrumental threshold of the peak flux (Pth) for a given 
burst with luminosity L 7 according to equation (11). The L J ul i n and L 7imax in our analysis 
are taken as 10 46 and 10 52 erg s _1 , respectively. 



3.2. Reference Sample for Constraining Model Parameters 

Free parameters in our models are po, L c (and its deviation <7£ ), 9 C (and its deviation 
cfq) and e. We take the current Swift/BAT GRB sample as a reference sample to constrain 
these parameters. Swift/BAT has detected 670 GRBs by June 2011. Among them 170 GRBs 
have redshift measurement. We bootstrap a mock Swift GRB sample and compare it with 
the observations. The peak fluxes, redshift and the spectral information for these GRBs are 
taken from published BAT catalogs (Sakamoto et al. 2008; Sakamoto et al. 2009; Sakamoto 
et al. 2011). The peak fluxes are measured in 1-s timescale. Note that the low- luminosity 
GRBs (L < 10 49 erg s -1 ) are not included since they may belong to a distinct population 
(Soderberg et al. 2004; Liang et al. 2007; Chapman et al. 2007). 



3.3. Bootstrapping Procedure 

(1) Generate a mock GRB characterized by z, L 7 , and 9y These parameters are ran- 
domly from the probability distributions of Eqs. (6), (9), and (10) for a given set of model 
parameters (L c , a Li , 9 C , <tq, e). Since z < 10 for the current BAT sample, the range of z for 
our analysis is from to 10. 

(2) Calculate the isotropic gamma-ray luminosity with L = L 7 /(l — cos^) and derive 
E' p in the burst frame with an L - E' p relation, i.e. E' p = 200 keV(L/10 52 ) 1 / 2 /C, where C is 
randomly distributed in [0.1, 1] (Liang et al. 2004). 

(3) Calculate the peak flux in the Swift/BAT band with Eq. 11 assuming that the photon 
indices in the energy bands lower and higher E p are -1 and -2.25, respectively (Preece et 
al. 2000), and compare it with the threshold of Swift/BAT. A mock GRB is recognized as 
detectable if its peak flux is over the threshold. 
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(4) Calculate the detection probability of a detectable GRB with the probabilities de- 
scribed by Equations (12), (13) and (14). 

(5) Repeat the above steps to make a mock redshift-known BAT GRB sample (170 
GRBs). 

(6) Evaluate the consistency between the mock and the observed sample by comparing 
the z, L, and log iV — log P distributions. We measure the consistency with the K-S test. We 
calculate the probabilities of the K-S tests, P£_s, Pr-s an< ^ Pr-s > f° r the l°g^ — logP, 
L and z distributions and then define the global K — S test probability as Pr_ s — Pr-s x 
Pr-s x Pr-s- A larger value of Pr-s indicates a better consistency. A value of Pr-s > 0.1 
is generally acceptable to claim the statistical consistency, whereas a value of Pr-s < 10~ 4 
convincingly rejects the hypothesis of the consistency. 

3.4. Results 

Following the procedure described above we search for the model parameters that can 
present the best consistency between the mock and the observed BAT GRB samples. In 
order to reduce the amount of calculation, we first randomly generate a large sample of the 
parameter sets(L c , ur , 9 C , <7g, e) in wide parameter ranges. Please note that we do not need 
to consider po m our bootstrap analysis since we make the comparisons between our results 
and observations with the relative probability distributions. We find that (<7l 7 ,<70, e) = (0.4, 
0.6, 0.4) can well reproduce the scattering of the GRB distributions in the L — (l + z) plane. 
We therefore fix the three parameters and then refine our search for the parameters L c and 
9 C . Figure (3a) shows the contours of the relative Pr_ s in the log L c — log 9 C plane. From 
Figure (3a) one can find that both 9 C and L c seem to be degenerate. If 9 C oc our results 
roughly reproduce the observations. This is reasonable based on the equation (11). It shows 
that for a given instrumental threshold P t h, a lower-L 7 GRB may be detectable if its 9- } is 
smaller. Therefore, one cannot constrain the intrinsic L c and 9 C based on a flux-limit sample. 
The parameter set (logL c /erg s _1 , log# c /rad) = (49.69,-1.27) gives the best consistency 
to the current Swift/BAT sample. Figure 4 displays the comparisons between our mock 
sample with the BAT sample in the logL — log(l + z) plane along with the distributions 
of log A" — logP, L, and z. The Prs values for the logL, log(l + z) and logA^ — logP 
distributions are also marked in each panel. One can observe that the observed BAT sample 
can be well reproduced with our bootstrapping method. In order to get a robust results, we 
also generate a large sample of 10 4 GRBs with the best parameter set and make the contours 
of the relative probability distribution of the GRBs in the log L — log(l + z) plane. The result 
is also shown in Figure 4, confidently suggesting a good consistency between our our results 
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and observations. 

The 9j— known sample has 77 GRBs. We randomly pick up a sub-sample of 77 GRBs 
from our mock GRB sample and show the 1- and 2-dimensional \og9j — log(l + z) GRB 
distributions in Figure 5. The K-S test probabilities are also marked in each panel in the 
figure. The K-S test yields a global K-S probability larger than 0.1, indicating a good 
consistency between the mock and observed BAT samples. As described in §3, there is no 
any intrinsical relation between Q } and z in our model. The observed Oj — z dependence in the 
mock sample definitely is due to instrumental selection effects. Therefore, the apparent 9 S — z 
dependence observed in the current GRB sample would be explained with the instrumental 
selection effect. 

As shown in Figure (3a), both 8 C and L c are degenerate. The triggered GRB number 
with an instrument in a given operation period may place constraint on the 6 C and L c , 
if the local GRB rate is known. However, the p is quite uncertain. It is reported that 
p ~ 1 Gpc^yr- 1 (Schmidt 2001; Lloyd-Ronning et al. 2004; Guetta et al. 2005; Liang et al. 
2007; Wanderman & Piran 2010). Please note that this rate is obtained without considering 
the jet collimation. Assuming a typical 9 } of 0.1 rad, the p is Po ~ 180 Gpc~ 3 yr~ 1 . The 
local Ib/c SN rate, which is reported as 4.8 x 10 4 Gpc _3 yr _1 (Marzke et al. 1998; Cappellaro 
et al. 1999; Folkes et al. 1999), presents a robust upper limit to p . BAT triggered 670 
GRBs in during 6.5 operation years. We therefore calculate po with equation 15 using the 
best parameter sets derived above. We get po = 2.85 x 10 2 Gpc~ 3 yr _1 for the parameter 
set (logL c /erg s _1 , log# c /rad) = (49.69,-1.27). The derived local GRB rates are higher 
than that reported in literature, but still under the local Ib/c SN rate. Our results suggest 
that ~ 0.59% of Type Ib/c SNe may be accompanied by a GRB. This may has profound 
implications for understanding the relation between GRBs and Type Ib/c SNe (e.g., Lamb 
et al. 2005). 



4. Conclusions and Discussion 

Assuming that the intrinsic distributions of both L 1 and 9j are Gaussian, we have 
shown that the observed 6*j dependence on redshift can be interpreted with the instrumental 
selection effects. 

The L 7 (or E 7 ) and 6>j are two basic characteristics of the GRB jets. As shown in Figure 
1, the observed distributions of 6- : and E 1 are log-normal. Frail et al. (2001) proposed the 
idea of a standard energy reservoir in GRB jets based on the narrow distribution of E 1 (see 
also Berger et al. 2003). Although our E 1 distribution is still log- normal, the width of the 
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distribution is much larger than that reported by Frail et al. (2001). Physically, it is difficult 
to access the intrinsic distributions of both L 7 and 8j from observations since the observed 
GRB sample with redshift measurement is heavily suffered from instrumental selection effects 
and observational biases (e.g., Bloom 2003). It was suggested that the intrinsic distributions 
of L 7 and #j are a power-law or a broken power-law (e.g., Lin et al. 2004). As shown by 
Qin et al. (2010), in case of a smooth broken power-law of isotropic peak luminosity L 
bootstrap analysis can well reproduce the observed L — z distributions in both the one- 
and two-dimensional planes. Intuitively, being due to both instrumental flux-limit and low 
triggered probability for low L 7 , the trigger probability for GRBs at low L 1 end would 
be reduced significantly. Similarly, the probability of a narrower jet pointing to the light 
of sight would be significantly lowered. Hence its detection probability would rapidly drop. 
Therefore, observationally, one cannot distinguish a Gaussian or a power-law distribution for 
with current GRB samples. In order to examine this expectation, we take the distributions 
of L 7 and 9 S as a power-law function, e.g., $(£ 7 ) = (L-,/L? L ) KL ; ip(9j) = (0j/9^ L ) Ke . We 
fix = 10 50 erg s _1 and 6^ L = 0.1 rad, but let kl and kq as free parameters. Our 
bootstrapping method still can well reproduce the observations. The contours of the relative 
Pks m the k>l — K e plane are also plotted in Figure (3b) for comparison with the P^s 
distribution in the 6 C — L c plane. Similarly, k l and n e seem to be degenerated as that of 6 C 
and L c . In case of Kg oc k°' 5 , our results roughly reproduce the observations. The parameter 
set (KL,Kg) = (—1.02,-1.42) gives the best consistency between our results and the BAT 
GRB sample, which is shown in Figs. 6 and 7. Yonetoku et al. (2005) suggested a k of 
~ —2 from CGRO/BATSE data, which is steeper than our result for the BAT sample. The 
derived p i n the parameter set (k Li kq) = (—1.02, —1.42) is 1.52 x 10 2 Gpc _3 yr _1 . 

Some apparent correlations between some GRB properties and redshift were discussed 
by some authors with small samples, such as the redshift dependence of the spectral lag of 
GRBs in prompt phases (Yi et al. 2008) and the duration of the shallow decay phase of GRB 
X-ray afterglows (Stratta et al. 2009). Physically, these correlations may be results from 
the cosmic evolution of the explosion energy and/or the jet geometry (e.g., Xu et al. 2005), 
if they are intrinsic. However, as shown in our analysis, there are no convincing cosmic 
evolution features. Therefore, these apparent correlations may be, at least partially, due to 
sample selection effects or observational biases. 

Tight correlations between E 1 and E' p (Ghirlanda relation; Ghirlanda et al. 2004) 
or among £i so , t- s and E' (Liang-Zhang relation, Liang & Zhang 2005) were discovered with 
small GRB sample. As shown in this analysis, the E 1 distribution would suffer great selection 
effects. It is known that the observed distribution of E p is also due to the selection effect. 
The observed narrow distribution of E p by CGRO/BATSE is significantly broadened by the 
observations with HETE-2 and BAT. Our analysis cannot figure out if the selections effects 
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on both Ej and E p distributions can be canceled out and evaluates whether these relations 
are intrinsic. 
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Fig. 1. — Distributions of the jet opening angle (panel a) and the geometrically-corrected 
gamma-ray energy (panel b) for the 77 GRBs in Table 1. The solid lines are the best fits 
to the data with a Gaussian function. 
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Fig. 2. — Jet opening angles as a function of redshift for the 77 GRBs in Table 1. The GRBs 
with tj measuring at different bands are marked with different symbols as indicating in the 
legend. The solid line is the best fit to all data. 
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Fig. 3. — Contours of the global K-S test probability in the L c — 9 C plane (panel a) and 
kl — kq plane (panel b). The dash lines denote 9 C oc L®- 5 and kq oc respectively. 
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Fig. 4. — Comparisons of two-dimensional logL-log(l + z) distributions and one-dimensional 
log L, log (1 + z) and log P (accumulative logiV-logP) distributions between the 170 Swift 
GRBs (solid dots and solid lines) and our results (open dots and dashed lines). The contours 
in the logL-log (1 + z) plane show the normalized probability distributions of the mock GRBs 
with different color scale lines. One dimensional K-S test probabilities for the comparisons 
are marked in each panel. 
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Fig. 5. — Comparisons of two-dimensional log0j-log(l + z) distributions and one- dimensional 
log 9j and log (1 + z) distributions between the 77 GRBs with opening angles (solid lines) 
and our results (open dots and dashed lines). The contours in the log^j-log (1 + z) plane 
show the normalized probability distributions of the mock GRBs with different color scale 
lines. One dimensional K-S test probabilities for the comparisons are marked in each panel. 
Other symbols in the log^j-log (1 + z) plane are the same as in Fig. (2) 
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Fig. 6. — The same as Fig 4 but for power law distributions of the jet-corrected luminosity 
and jet opening angle. 
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Fig. 7. — The same as Fig 5 but for the power law distributions of the jet-corrected luminosity 
and jet opening angle. 
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060605 


3 


8 


0.70±0.09 


15-150 


43.6 


0.24±0.02 


o 





027+0.001 






44,44,57 


060614 





13 


28.20±0.40 


15-350 


2.04 


1 45+ 1 - 16 
4-^° — 0.41 


X 





,,0 + 0.040 
— 0.014 


{1 


018 


44,46,45 


060707 


3 


43 


1.60±0.15 


15-150 


79.9 


12.26+\ 2 7 \ 


X 





,,,+0.018 
J-J-J- q Q20 





495 


44,44,45 


060714 


2 


71 


2.83±0.17 


15-150 


109 


0.12±0.01 


X 





020+0.001 





022 


44,44,58 


060729 





54 


2.61±0.21 


15-150 


4.70 


86.23+g" 1 


X 





314 +0.155 
° -0.028 





229 


44,44,45 


060814 





84 


14.60±0.24 


15-150 


88.7 


0.55±0.14 


X 





048+0.005 





101 


44,44,45 


060906 


3 


69 


2.21±0.14 


15-150 


138 


16+ 004 
uiD -0.03 


X 





non + 0.002 
u ^ u -0.001 





028 


44,44,45 


060908 


2 


43 


2.80±0.11 


15-150 


101 


01+ 02 
u,ul -0.003 


X 





nno + 0.005 
UU8 -0.001 





004 


44,44,45 


060926 


3 


21 


0.22±0.03 


15-150 


11.5 


0.06±0.05 


X 





019+0.006 





002 


44,44,45 


060927 


5 


6 


1.13±0.07 


15-150 


119 


0.04±0.02 


X 





011+0.002 





007 


44,44,45 


061121 


1 


31 


13.70±0.20 


15-150 


280 


34+ ' 18 


X 





nqo + 006 

"■^^-o.oos 





140 


44,44,45 


070125 


1 


547 


22.50±3.50 


20-10000 


108 


1 ()5 + 035 
1 ' uo -0.30 


X 





()r „ + 0.007 
uoo -0.006 





149 


59,59,45 


070208 


1 


17 


0.45±0.10 


15-150 


3.66 


11+ 005 
u,11 -0.03 


X 





no 7 +0.007 
UJ '-0.004 





003 


44,44,45 


070306 


1 


5 


9 ()0 +0 - so 
J ' uu -0.80 


15-350 


88.8 


, ,, + 1.79 
iJO -086 


X 





ngg + 0.030 
uoa -0.014 





156 


44,46,45 


070318 





84 


2.48±0.11 


15-150 


13.2 


, -_+2.88 
°-° -0.63 


X 





,,, + 0.037 
lz ^-0.008 





098 


44,44,45 


070411 


2 


95 


2.70±0.16 


15-150 


115 


24+°' 6s 


X 





„o 5 + 0.023 
U ' !;i -0.005 





037 


44,44,45 


070508 





82 


19.60±0.27 


15-150 


124 


58+ 93 
u - oo -0 .50 


X 





„ 47 +0.028 
-0.015 





136 


44,44,45 


070611 


2 


04 


0.39±0.06 


15-150 


9.25 


1 13 + 0-39 


X 





069 +0.009 
UD -0.012 





022 


44,44,45 


070714B 





92 


0.72±0.09 


15-150 


19.2 


0.01±0.002 


X 





013+0.001 





002 


44,44,45 



-25 - 



Table 1 — Continued 



GRB 


z 




Energy" 


E iso (~f) b 


H 


Band c 


ef 


E -< 


Refs. 






(10~ 6 erg cm~ 2 ) 


(keV) 


(10 61 crg) 


(days) 




(rad) 


(10 51 erg) 




070721B 


3.63 


7.40±0.60 


15-350 


365 


n + 001 
u,±± -0.02 


X 


0.015±0.001 


0.043 


60,60,45 


070810A 


2.17 


0.69±0.06 


15-150 


17.3 


09+ 04 
u.ua_ 06 


X 


024+ 004 


0.005 


44,44,45 


071003 


1.1 


8.30±0.30 


15-150 


206 




X 


037+ 002 
U UJ ' -0.003 


0.139 


44,44,45 


071010A 


0.98 


0.20±0.04 


15-150 


1.22 


„ „,+0.22 
0.81_ Q 2Q 


X 


092+ 010 
U.UM2_ oog 


0.005 


44,44,45 


071010B 


0.947 


4.40±0.10 


15-150 


19.1 


3.44±0.39 


O 


0.113±0.005 


0.121 


44,44,61 


071031 


2.69 


0.90±0.13 


15-150 


34.3 


7l+ - 35 
u - ' 1 -0.50 


X 


046+ 008 

U.UTO_ Q 012 


0.036 


44,44,45 


080319B 


0.937 


81.00±1.00 


15-150 


1173 


0.03±0.01 


O 


0.012±0.002 


0.080 


44,44,62 


090323 


3.568 


202 00+ 28 00 


20-10000 


3722 


17 80+ 19 60 
i I .nu_ 1 70 


X 


078+ 032 
u - u ' o -0.003 


11.4 


63,63,64 


090328 


0.736 


95.00±10.00 


8-40000 


100 


6 40+ 12 00 


X 


n ,,,,+0.085 


0.729 


65,65,64 


090902B 


1.822 


374.00±3.00 


50-10000 


2972 


a ~n + 2.-10 


o/x 


n ncc + 0.009 

0.065l Q 003 


6.27 


66,67,64 


090926A 


2.106 


145.00±4.00 


8-1000 


1808 


9.00±2.00 


o/x 


0.077±0.006 


5.32 


68,69,64 



a The energy range (keV) over which the fluence was reported. 
b The isotropic gamma-ray energies (Ej so (~f)) 

c Thc jet-break is identified in radio (R), optical (O), or X-ray (X), respectively. 
^The jet opening angle (#j) 

e The geometrically-corrected gamma-ray energy (E~ f ) 
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